#include "Multigrid.h"
using namespace std;

int main(){
    int n = 128,M = 20;
    //M:the number of maximum iteration
    double u0_[n+1],un_[n+1],u_0[n+1],u_n[n+1],epsilon;
    //u0,un: boundary conditions

    double h = 1.0/n;
    double v0[(n-1)*(n-1)];
    double f0[(n-1)*(n-1)];
    double u[(n-1)*(n-1)];
    u00 = 1;
    u0n = 1;
    un0 = 1;
    unn = exp(sin(1)*sin(1));
    for(int i = 0;i < n-1;i = i + 1){
        double x = (i+1)*h;
        for(int j = 0;j < n-1;j = j + 1){
            double y = (j+1)*h;
            v0[i*(n-1)+j] = 0;
            f0[i*(n-1)+j] = exp(sin(x)*sin(y))*(-2*sin(x)*sin(y)+cos(x)*cos(x)*sin(y)*sin(y)+cos(y)*cos(y)*sin(x)*sin(x));
            u[j] = exp(sin(x)*sin(y));
        }
    }
    for(int i = 0;i < n+1;i = i + 1){
        double x =
        u0_[
    }
    f0[0] = f0[0] + u0/h/h;
    f0[n-2] = f0[n-2] + un/h/h;
}
